Striped phases in the two-dimensional Hubbard model with long-range Coulomb 

interaction 
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We investigate the formation of partially filled domain walls in the two-dimensional Hubbard 
model in the presence of long-range interaction. Using an unrestricted Gutzwiller variational ap- 
proach we show that: i) the strong local interaction favors charge segregation in stripe domain 
walls; ii) The long-range interaction favors the formation of half-filled vertical stripes with a period 
doubling due to the charge and a period quadrupling due to the spins along the wall. Our results 
show that, besides the underlying lattice structure, also the electronic interactions can contribute 
to determine the different domain wall textures in Nd doped copper oxides and nickel oxides. 



PACS numbers: 71. 27,71. 10.Fd,75.10.Lp,75.60.CH 

The occurence of charged domain walls in the high- 
T c superconductors presently attracts a lot of interest 
also with regard to possible pairing scenarios In- 
commensurate spin correlations have first been observed 
in La2- K Sr a: Cu04 (LSCO) by neutron scattering exper- 
iments |[|-[|- More recently it was found that the in- 
commensurate spin fluctuations are pinned in Nd doped 
LSCO and nickel oxide compounds leading to spin- 
and charge-stripe order in these materials p0| . However, 
whereas the domain walls in the hole-doped La2Ni04 sys- 
tem are oriented along the diagonals of the MO2 lat- 
tice, it turns out that in Lai.4sNdo.4Sro.i2Cu04 the ori- 
entation of the stripes is along the Cu-0 bond direction. 
Moreover, the hole concentration in the domain walls is 
one hole per Ni site in the nickelates and one hole ev- 
ery second Cu site in the Nd-doped cuprates. A com- 
parison of the low-temperature orthorhombic and the 
low-temperature tetragonal structure suggests that the 
vertical stripes in La1.48Ndo.4Sro.12Cu.O4 are pinned by 
the tetragonal lattice potential whereas the orthorombic 
phase in the nickel oxides favors a diagonal orientation. 
We will show here that the electronic interactions may 
also play an important role in establishing different do- 
main wall structures. 

The stripe instability for doped antiferromagnets was 
predicted theoretically in JllJ within Hartree-Fock (HF) 
theory applied to the extended Hubbard model and con- 
firmed by a number of subsequent investigations |fl2f . For 
small values of the Hubbard on-site repulsion U (gener- 
ally smaller than 3t - At) these calculations result in a 
striped phase oriented along the (10)- or (Ol)-direction 
whereas for higher values of U the orientation is along 
the diagonals. Within HF theory the stripe solutions 
become unstable for U > 8t towards the formation of 
isolated spin-polarons. 

However, all stripe calculations performed so far within 
the HF approximation of the Hubbard model predict one 
hole per site along the domain wall (filled stripe) p3|. 



This contrasts the observation of half-filled stripes (i.e. 
with half a hole per site) in the Lai.6- a ;Ndo.4Sr : j : Cu04 
system which is a yet unresolved problem of mean-field 
theory. In Ref. it was addressed the question whether 
the inclusion of an additional nearest neighbor repulsion 
V in the Hubbard-hamiltonian may favor the formation 
of partially filled stripes. According to Ref. |Q the half- 
filled stripe solution is stabilized by a quadrupling of the 
charge- or spin-period along the stripe. However, al- 
though the nearest neighbor repulsion slightly enhances 
the stability of the half-filled wall this never corresponds 
to the HF ground state for realistic parameter values. 
Instead the main effect of V is to shift the crossover to 
isolated spin polarons to lower values of U. 

In the present paper we show that a proper treat- 
ment of the strong local repulsion U plays an indirect 
but crucial role in stabilizing half-filled vertical domain 
walls. Specifically we apply a slave-boson version of the 
Gutzwiller approach within an unrestricted variational 
scheme. Contrary to the pure HF approach, which heav- 
ily underestimates the effective attraction between the 
charge carriers and predicts repulsion for U > 6t, it was 
recently shown that within the slave-boson scheme the 
attraction persists up to very large U As a con- 

sequence of this more suitable treatment of the strong 
coupling limit, U greatly favors the charge segregation 
in striped domains as opposed to spin polarons. In the 
absence of long-range (LR) forces, completely filled di- 
agonal stripes stay more stable than half-filled vertical 
ones. However, due to their increased stability with re- 
spect to isolated polarons, the stripe solutions now allow 
for a less disruptive introduction of stronger LR forces, 
which affect the completely filled stripes more than the 
half- filled ones. Then, for a sizable but still realistic LR 
repulsion, the half-filled vertical stripe may become the 
ground-state configuration. 

We consider the two-dimensional Hubbard model on a 
square lattice, with hopping restricted to nearest neigh- 
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bors (indicated by the bracket < i,j >) and an additional 
LR interaction: 

H = -t 4<r c i,cr + U'Y2 n iA n hl + E v ii n i,<? n i,<?' 

(1) 

where destroys an electron with spin a at site i, and 
ni,a = c\ a Ci. a - U is the on-site Hubbard repulsion and 
t the transfer parameter. For all calculations we take 
t = 1. For the Coulomb potential we assume an interac- 
tion of the form Va — , = where the param- 

eters Vq and a are specified through the on-site repul- 
sion U and the nearest-neighbor interaction V^ i?t +i via 

a = V 7hn+ i/^Ju 2 - V 2 n+1 and V = all/2. Since we 
consider a finite lattice with periodic boundary condi- 
tions it is also necessary to restrict the LR part to half 
of the lattice dimension both in x- andy-direction. 

Following Kotliar and Ruckenstein [|16| we enlarge the 
original Hilbert space by introducing four subsidiary bo- 
son fields ei, Sirf, s^x, and di for each site i. These op- 
erators stand for the annihilation of empty, singly oc- 
cupied states with spin up or down, and doubly oc- 
cupied sites, respectively. Since there are only four 
possible states per site, these boson projection opera- 
tors must satisfy the completeness constraints eje^ + 

+ d\di = 1 and n i:U = s\,a s i,° + A^- In 
the saddle-point approximation, all bosonic operators are 
treated as numbers. Furthermore we will approximate 
the LR part by a HF decoupling. The resulting effec- 
tive one-particle Hamiltonian describes the dynamics of 
particles where the hopping amplitude between states 
(i,er) and (j,c) is renormalized by a factor z\ G Zj.o with 

z%,o = (ef + sl_ a y 1/2 (ejS i)(T + Si- a di) (d 2 + s 2 i a y 1/2 
The effective one-particle Hamiltonian can be diagonal- 
ized by the transformation a, a = ^2 k &i.a{k)cik where 
the fermion wave functions <E>i iCr (&;) obey the orthonor- 
mality constraint £\ ff = <W 

Given a system with N e i particles we obtain for the 
total energy 

— , Nal — , 

E tot = -t ]r <^>E $ ^( fc )^>( fc ) +C7 E d " 

<ij>,a k—1 i 

+ E v ij ^*l a (k)*iAkWAk')*i,Ak') (2) 

i^j,a(T ! k,k' 

which has to be minimized with respect to the fermionic 
wavefunctions and bosonic fields within the orthonormal- 
ity and completeness constraints (for further details of 
this approach see fq]). 

Single-stripe calculation. - In the presence of LR in- 
teractions a single stripe will always result unstable with 
respect to isolated polarons by increasing the length of 
the stripe. Indeed the Coulomb energy per hole of a 



charged wall of length L increases as log L and it is not 
compensated by the coupling to a uniform distribution 
of background charges of opposite sign, as it is the case 
for a regular array of stripes. Nevertheless the analysis 
of the single-stripe case allows to illustrate the stabiliza- 
tion effect of a large U (within the Gutzwiller approach 
and with respect to the HF treatment) and to find out 
the most stable stripe configuration within each class of 
stripes. 

It has been intensively discussed in Ref. Il4j that for 
the half-filled wall to become a saddle-point for any effec- 
tive one-particle description, it is necessary to quadruple 
the period along the stripe. In fact there exist several 
possibilities of performing this period quadrupling and 
in Fig. 1 we sketch the charge- and spin structures of the 
two types of half-filled vertical domain walls which have 
the lowest energy among the different realizations. We 
find that the stripe configurations investigated in Ref. 
|l4| (cf. Figs. 8,9 therein) are always higher in energy 
than the structures shown in Fig. 1. 
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FIG. 1. Charge- and spin-density structure for the two 
half- filled vertical domain wall types described in the text. 
The diameter of the circles indicates the hole charge and the 
length of the spin arrows symbolizes the value of (S z ). 

Fig. 1(a) is of the CDW type where the main charge 
modulation is with q?, = n along the wall whereas the 
spin varies with qf, = tt/2. We find that among the half- 
filled vertical walls this configuration is lowest in energy 
in the range At < U < 7.5t. However, for strong on-site 
repulsion (U > 7.5t) the charge and spin realization with 
the lowest energy corresponds to the staggered structure 
indicated in Fig. 1(b) although the difference in energy to 
the CDW- type of Fig. 1(a) is rather small (the difference 
in energy per hole for U — lOi is w 5 * 1CP 3 <). 

For the single stripe solutions the size of the supercell is 
9x8. Calculation of the diagonal structures for periodic 
boundary conditions without frustrating the background 
spin configuration in principle requires an 'uneven x un- 
even' lattice which does not allow for the investigation 
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of half-filled diagonal walls. For this reason and to avoid 
the comparison between different lattice sizes for differ- 
ent domain wall types we have choosen supercells with 
periodic boundaries for vertical stripes and with shifted 
periodic boundaries for diagonal stripes. In the 'shifted 
boundaries' the equivalent supercells in the vertical di- 
rection are horizontally shifted by one lattice spacing (cf. 
Fig. 4 in Ref. p^|). To avoid double counting in the 
Coulomb interaction energies, we cut off the LR forces at 
"half-minus-one" the size of the supercell. 

The energy per hole of each configuration was calcu- 
lated with respect to the reference state of a uniform 
AF lattice with one particle per site: E(config) = 
(E(N h ) - E AF (0)) /N h 0. In particular the energy 
of an isolated spin-polaron is given by E{polaron) = 
E(Nh = 1) — Eaf(0)- To compare the stability of stripes 
on our finite-size system, we report the energy differences 
E(config) — E(polaron). 

The results are displayed in Fig. 2 as a function of the 
nearest neighbor value of the Coulomb repulsion V, lin +i 
for U = 7t in Fig. 2(a) and for U = lOt in Fig. 2(b). In 
case of half-filled vertical stripes the curves in Fig. 2(a) 
and 2(b) correspond to the structures in Fig. 1(a) and 
1(b) respectively. 




b) 








/ srf 

V / 
/ p 


V 


/ / 


U=10t 



0.0 0.1 0.2 0.3 0.4 
V V 

FIG. 2. Energy difference of various striped phases and the 
isolated polaron lattice E(stripe) — E(polaron) as a function 
of the nearest neighbor contribution V n , n +i to the long-range 
potential. Solid line (O): completely filled vertical stripe; 
long-dashed line (+): half-filled vertical stripe; short-dashed 
fine (□): completely filled diagonal stripe; dotted line (x): 
half-filled diagonal stripe. System size: 9 x 8; a) U=7t; b) 
U=10t. 

Disregarding the eventual above-mentioned instability 
of the single stripes with respect to isolated polarons 
by increasing lenght, various features are worth not- 
ing. First of all, completely filled stripes increase more 
rapidly their energy (have a larger slope) than half-filled 
ones upon increasing the LR repulsion. Therefore LR 
forces favor half-filled stripes, which eventually become 
the most favorable wall textures. 

However, the most relevant effect to be noticed here is 



the role of a large local repulsion U affecting the ener- 
gies of the various textures. In particular a comparison 
between Fig. 2(a) and Fig. 2(b) shows that U strongly 
reduces the energy of the half-filled stripes with respect 
to the filled ones already at V — 0. By combining this 
reduction with the effects of LR forces on the stripes 
it follows that increasing U makes the half-filled vertical 
stripe the most stable among the wall solutions at smaller 
values of V n , n +i- 

Interstripe interaction and spin-polaron lattice. - To 
assess the actual ground state we now consider the inter- 
action between stripes or between spin polarons. Indeed 
one cannot neglect the interstripe repulsion since, for ex- 
ample, in Lai.4gNdo.4Sro.i2Cu04 the stripe separation is 
four times the Cu-Cu distance only. To incorporate the 
repulsion between stripes, we have calculated the energy 
of vertically oriented stripes on a 32 x 4 lattice. For a 
concentration of 1/8 this results in an array of 4 com- 
pletely filled or 8 half-filled stripes. The energy of these 
arrays with respect to the polaronic Wigner lattice are 
depicted in Fig. 3 (the Wigner lattice now corresponds 
to 16 spin polarons with maximum distance). 




FIG. 3. Energy difference between vertically striped phases 
and the 'interacting' spin-polaron Wigner lattice as a function 
of the nearest neighbor contribution V n , n +i to the long-range 
potential (in units of t). Solid line (O): completely filled ver- 
tical stripe; long-dashed line (+): half-filled vertical stripe. 
System size 32 x 4, U=10t. 

In this case we obtain a crossover to half-filled vertical 
domain walls for V ~ 0.12i. It is interesting to observe 
that now the stripe solutions gain in energy upon switch- 
ing on the LR part with respect to the polaron lattice. 
This enhances the parameter range of stability for the 
stripes, which no longer become unstable towards the 
decay into isolated polarons. 

To assess the absolute stability of the half-filled verti- 
cal stripes, we should also compare their energy with the 
filled diagonal stripes, which, in the single-stripe analysis 
result to be more stable than the filled vertical stripes. 
However, filled diagonal stripes are strongly destabilized 
by the elongated shape of the 32 x 4 supercell so that 
we do not included their energy in Fig. 3. As an alter- 
native to the direct calculations on the elongated clus- 
ter, to extract informations about diagonal stripe config- 
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urations, we analyzed a two-dimensional regular array of 
charged wires with fixed global charge density. We have 
found that the electrostatic potential energy is lower for 
wires with higher linear charge density at a larger dis- 
tance than for less charged wires more closely spaced. 
Therefore diagonal stripes, which at given planar den- 
sity are closer by y/2, but less densely charged by the 
same factor, are less favorable than the vertical stripes 
as far as the electrostatic Coulombic energy is concerned. 
On the other hand, our single-stripe investigation already 
demonstrated that a proper treatment of the strong local 
repulsion U opens the way to a stabilization of half-filled 
vertical stripes with respect to the filled diagonal stripes. 
From the above purely electrostatic analysis and from 
the results of Fig. 3 we can therefore safely conclude 
that half-filled vertical stripes are the ground-state con- 
figuration for V > 0.12i at large enough U. 

To summarize, we have shown that a LR Coulomb in- 
teraction added to the 2D Hubbard model gives rise to 
half-filled vertically oriented domain walls when treated 
within an unrestricted Gutzwiller approach. This fea- 
ture does not appear in semiclassical Hartree-Fock ap- 
proximations where the effective attraction between spin- 
polarons is underestimated. Depending on the value of 
the on-site repulsion U, we expect the domain wall struc- 
ture in the Nd-doped LaCuO system to be of the type 
shown in Fig. 1(a) and Fig. 1(b), respectively. More- 
over, the here investigated competition between com- 
pletely filled diagonal and half-filled vertical stripes can 
explain the different hole orderings in nickelates and in 
Nd cuprates even without invoking a relevant role of lat- 
tice interactions. Specifically, our findings suggest that 
nickelates could be characterized by a smaller U < 7t 
and/or a smaller V accounting for their filled diagonal 
stripes. On the other hand, although the half-filled ver- 
tical stripes in the Nd-enriched LaSrCuO systems are to 
some extent likely fixed along the (1,0) direction by the 
underlying lattice structure, we showed here that, if large 
values of U and sizable values of V characterize these sys- 
tems, then electronic correlations would also contribute 
to give rise to vertically half-filled stripes. 

In general the filling and orientation of the stripes de- 
pend on the specificity of the electronic forces and struc- 
tures and therefore inside the various oxide families dif- 
ferent textures of the stripe phase may prevail. 
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